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Small RNA(sRNA)-mediated post-transcriptional regulation differs from protein-mediated regulation. Through base- 
pairing, sRNA can regulate the target mRNA in a catalytic or stoichiometric manner. Some theoretical models were built for 
comparison of the protein-mediated and sRNA-mediated modes in the steady-state behaviors and noise properties. Many 
experiments demonstrated that a single sRNA can regulate several mRNAs, which causes crosstalk between the targets. 
Here, we focus on some models in which two target mRNAs are silenced by the same sRNA to discuss their crosstalk 
features. Additionally, the sequence-function relationship of sRNA and its role in the kinetic process of base-pairing have 


been highlighted in model building. 
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1. Introduction 


In bacteria, small RNAs (sRNA) post-transcriptionally 
regulate target mRNAs through base-pairing, which plays 
an important role in physiological processes, such as iron 
metabolism,!!] acid stress, 2] quorum sensing, 1 and outer 
membrane protein expression. 4! sRNA have been categorized 
into two types according to the location in the genome: cis- 
encoded sRNA and trans-encoded sRNA. The former is lo- 
cated in the same region as its target mRNA in the genome 
and completely bases pairing with the target, while the lat- 
ter is located in regions different from those of its targets and 
partially bases pairs with the targets.! Here, we focused on 
trans-encoded sRNAs, which usually have three basic regions: 
the core interaction region, the Hfq-binding site, and a rho- 
independent terminal. The core interaction region, also called 
the “seed ” region, helps sRNA to specifically discern target 
mRNAs and initiates the base-pairing with mRNAs.'°! The 
Hfq-binding site usually an AU-rich region, helps sRNA bind 
at the proximal face of Hfq.!’! Hfq is a donut-shaped homo- 
hexamer that facilitates the base-pairing of sRNA and mRNA 
through improving sRNA stabilization or structural transfor- 
mation of mRNA.|®:°] However, the exact mechanism under- 
lying the effect of Hfq is not clear. Recently, some exper- 
iments have revealed that polyU of the terminal of sRNA 
is another Hfq-binding site and is necessary for the efficacy 
of sRNA.!0!1 SRNA has two effects on the translation of 
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mRNA: negative regulation, in which sRNA base pairing with 
mRNA near the ribosome binding site (RBS) inhibits the 
loading of ribosomes, and thus the translation of mRNA is 
repressed, [12,13] and active regulation, [14-16] in which the bind- 
ing of sRNA loosens the intra-base-pairing of RBS of mRNA 
so that the ribosome can carry out the translation. In the sce- 
nario of negative regulation, most of the sRNA-mRNA duplex 
is stoichiometrically degraded by endonuclease.!!”! Various 
models of sRNA-mediated regulation are shown in Fig. 1. 

Recently, another model has demonstrated that Chix 
(MicM) sRNA catalytically regulates chip (ybf{M) mRNA, 
which indicates that sRNA does not undergo codegrada- 
tion with mRNA in a stoichiometric manner but with 
recycling.!!*-!°l Based on the results of these experiments, 
Hao et al. built a catalytic-sRNA-mediated gene silencing 
(CSMGS) model to analyze its novel properties. 7°! For trans- 
encoded sRNAs, they usually regulated several mRNAs that 
comprise a network, [21,22] such as ChiX sRNA which is also 
codegraded with another target chbB mRNA, and is exhausted 
in a stoichiometric model to release chiP mRNA.!!%?3! 

With the accumulation of experimental data, several mod- 
els have been built to characterize the regulatory system. In 
this review, first, we have presented an sRNA-mediated reg- 


r,!74] and compared 


ulatory model in a stoichiometric manne 
it with the protein-mediated post-transcription model. Sec- 
ond, we have discussed the features of the catalytic model and 


stoichiometric model with respect to sRNA-mediated regula- 
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tion from both steady-state analysis and noise analysis. Third, 
we have introduced a model depicting the simplest network 
of two target mRNAs regulated by a single sRNA, which 
shows the hierarchical crosstalk behavior between the target 


mRNAs. 


the sequence-function relationship of sRNA in gene silencing, 


Fourth, we have discussed the latest research on 


which should provide the foundation and basis of modeling 
the sRNA-mediated regulation. 
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Fig. 1. (color online) Small RNA-based regulatory mechanisms of mRNA expression. (a) sRNA blocks translation of mRNA by directly 
binding to the ribosome binding site (RBS) of mRNA. (b) and (c) The binding of sRNA on mRNA sequesters the ribosome, because of which 
endonuclease sites are exposed to degradation by RNase E in a passive or active manner. (d) The RBS region of mRNA shows intra-base-pairing 
in the absence of sRNA that inhibits translation of mRNA. The binding of sRNA loosens the RBS so that it is available in a form of a single 
strand for translation by the ribosome. The grey ovals indicate the ribosomes that are involved in the translation of mRNA. The purple circles 
indicate the proteins generated. AUG is the start codon for translation of mMRNA.U7] 


2. Comparison of sRNA-mediated and protein- 
mediated regulation 


sRNAs post-transcriptionally regulated the expression of 
mRANs, which differs from the transcriptional regulation of 
transcription factors (TFs), a kind of protein-mediated reg- 
ulation. Levine et al. found that sRNA-mediated regula- 
tion exhibited novel properties, such as the threshold-linear 
property and that it more tightly repressed target mRNA be- 
low the threshold and showed better robust noise resistance 
than protein-mediated regulation.*! For further discussion 
of the characteristics, we introduce the work of Mehta et 
al.; they treated TF-mediated and sRNA-mediated regulation 
as a signal-processing system and accordingly discussed the 
steady-state behavior, noise properties, frequency-dependent 
gain, and dynamical response to the input signals. [5] 

Mehta et al. built a model for sRNA-mediated regulation 
by using mass action equations combined with intrinsic noise 
defined by Langevin terms as fellows: 


ds 1 
= Os Ty S 


ums t ts + Au, 


dm =] A A 

ad = An —T, M— UMS + n+ u, 

d = A 

+ = aym—t,|p + ñp. (1) 


In this model, the sRNA (mRNA) transcription rate is set 
as Œs (Œm), and degradation rate is set as a (Tn). L is the 
binding strength of sRNA and mRNA. The protein is trans- 
lated from mRNA at the rate œp. s, m, and p refer to the 
amount of sRNA, mRNA, and protein, respectively. fs, Am, 
fu, and ÎÌp characterize the noise related to the creation and 
degradation of sRNA, mRNA, sRNA-mRNA, and protein, re- 
spectively. The Langevin terms are formatted by the time cor- 
relation functions (j = s, m, p, U): 


MONE) = oF 5(t—-1'). 


The TF-mediated transcriptional regulation is modeled 


using the Langevin equations: [26-29] 
dm sdf 7 
dt = On — Ty m+n, 
dp 7 : 
y om Ty P+ Ap, (2) 
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where fm and f}, model noise in the synthesis and degrada- 
tion of mRNA and protein respectively, and obey the equations 


(j =m, p) 
(MONE) = 27; "FS —1'). 

Ignoring the Langevin terms and setting the derivatives 
to zero in Eqs. (1) and (2), we can determine the amount of 
proteins in the steady state in both sRNA-mediated and TF- 
mediated regulation. Consistent with the study of Levine et al., 
sRNA-mediated regulation shows the linear-threshold prop- 


erty for the target mRNA expression level with the threshold 
around Qn ~ d;.'24! When Œm < Os, the translation of mRNA 


is repressed; when Œm >> ds, the expression level of protein 
is linear to (Œm — Os). The binding strength u just determines 
the smoothness of the crossover around the threshold œs. By 
contrast, the amount of mRNA in the TF-regulated steady state 
is a linear function of Qn. [2029] 

Based on the equations above, they also studied the 
intrinsic noise properties,°°?! frequency-dependent gain 
(amplification), 2?! and dynamical response to large input sig- 
nals of sRNA-mediated and protein-mediated regulation. 4! 
Taken together, their characteristics are summarized in Table 


1 cited from Mehta et al. 251 


Table 1. Summary of the advantages and disadvantages of sRNAs as compared with transcriptional protein regulators (transcription factors). 


Advantages of sRNAs 


Advantages of protein regulators 


e sRNAs are better than protein regulators at keeping proteins “off” 
because a large pool of sRNA shortens the effective mRNA lifetime and 
buffers against target mRNA fluctuations. 

e sRNAs can filter high-frequency noise without compromising the 
ability to rapidly respond to large changes in input signals. 

= sRNAs likely fill a “niche” in allowing cells to transition quickly 
yet reliably between distinct states. 


e The intrinsic noise for sRNA based regulation is much larger than 
that for transcriptional regulation in a large (intermediate to high) range 
of expression levels of the regulated protein, especially in the presence 
of transcriptional bursting. 

e Protein regulators are better than sRNAs at transducing small input 
signals. 

=> Protein regulators are likely better suited for quantitative adjustment 
of protein expression than sRNAs. 


3. Models of sRNA-mediated gene silencing: 
Stoichiometric and catalytic modes 
3.1. Steady-state behavior of target mRNA regulated by 
sRNA 
sRNA-mediated gene silencing was initially believed to 
occur in the stoichiometric mode, which differs from the 
protein-mediated catalytic mode. However, some experi- 
ments in bacteria showed that sRNA can exert its effects in 
the catalytic mode also, which means that sRNA can pro- 
mote the degradation of mRNA with a certain ratio of sRNA 


survival. |1819] 
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Fig. 2. (color online) Catalytic-sRNA-mediated gene silencing 


(CSMGS) model. 2° 


To investigate the characteristics of sRNA-mediated reg- 
ulation with recycling, Hao et al. developed the catalytic- 
sRNA-mediated gene silencing model (CSMGS),!°! based on 
Levine’s model (Fig. 2).4! They introduced the recycling 
rate 6 referring to how much sRNA is surviving from the 


sRNA-mRNA combo degradation. The resulting recycling 


item OBs,(ssn) was added in the catalytic model shown in 
Eqs. (3), in which Psm(s-m) means the degradation rate of 
sRNA-mRNA combo. Therefore, setting 6 to 0 transforms 
the catalytic model into the stoichiometric model. As men- 
tioned above, Hfq is necessary for most of the trans-encoded 
sRNA regulation, but the mechanism is not very clear. Here 
both models assumed that Hfq is abundant. 55! 
The mass action equations of CSMGS model are 


- = A; — Bss —k-s-m+ bByn(sm), 
Stem) = k-s-m— Bsm(sm), 
z = kpm — Bpp. (3) 


In this model, s (m), Œs (Qn), and Ps (Bm) indicate the 
sRNA (mRNA) concentration, the transcription rate, and the 
degradation rate, respectively. The binding rate k refers to the 
strength of sRNA-mRNA interaction which is relevant to the 
kinetic process of base-pairing of sRNA with mRNA. 

Setting the differential equations to zero yields 


Os — Bss —k(1 — ô) -s-m=0, 
On — Bmm—k-s-m=0. (4) 


When 0 < 6 < 1, the stationary solution of sRNA, 
mRNA, sRNA-mRNA duplex, and protein can be obtained as 
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Fig. 3. Stationary solutions and Fano factors of target mRNA vary 
with its transcription rate Œm in different recycling ratio in the CSMGS 
model. 701 


The parameter À is defined as A = Pmps/k, which in- 
dicates the ratio of the spontaneous and codegradation rates. 
n = kp/Pp refers to the burst size of target protein produc- 
tion. Some experiments have suggested that the base-pairing 
of sRNA with mRNA is strong and fast; therefore, we can ap- 
proximate A to zero. 20! Then the stationary solution of mRNA 
can be written as follows: 


™=% ow, —a/(1—8) P (6) 


Figure 3(a) shows the steady-state expression level of 
mRNA responding to its transcription rate for different recy- 
cling ratios of sRNA. In the case of ô = 0, the CSMGS model 
refers to the stoichiometrical regulation mode of sRNA, whose 
threshold-linear property has been discussed in the model of 
Levine et al.?4 The CSMGS model exhibited the similar 
threshold-linear property, but with the threshold characterized 
by a;/(1 — 6). When the transcription rate of mRNA crosses 
the threshold, the amount of mRNA is linearly related to Qn. 
Compared to the stoichiometric mode, in the catalytic mode, 
gene silencing by sRNA is more rapid because a lesser number 
of sRNA are required to degrade mRNA with a given Œm. 206] 
3.2. Intrinsic noise analysis for the models of sRNA- 

mediated gene silencing 


It is believed that noise significantly affects the dynam- 
ics behavior of molecules in vivo such as sRNAs, mRNAs 


[37.38] There are mainly two types of noise: in- 


and proteins. 
trinsic noise, which is related to the low reactant number and 
reactant stochastic bursting, 5850] and output noise, which is 
related to diffusion between reactant pools and the neighbor- 
ing environment.'*"] Here we focus on the intrinsic noise in 
sRNA-mediated gene silencing. 

Generally, there are three approaches to study the noise of 
the biochemical reactions. The first approach is by means of 
Langevin equations (as in Eq. (1)), 44 which was further de- 
veloped as the linear noise approximation (LNA).&?3! This 
method can provide analytical results for both the mean value 
and variance of each reactant. The second approach is to use 


n.!38:441 All biochemical reactions can be 


the master equatio: 
depicted in the master equation, which gives the evolution of 
the joint probability of all the reactant species. Generally, it 
is difficult to solve the master equation analytically; therefore, 
many researchers will choose the third approach of stochastic 
simulation with the Gillespie algorithm. [454] 

One can choose the appropriate strategy according to 
the research objective; all these objectives can provide the 
mean value X and variances 0? of the reactant species (x = 
s,m, p,...), and the Fano factor Fy = o? /x, which signifies the 
noise strength of each biochemical reactant. In the CSMGS 
model, Hao et al. applied the Gillespie algorithm to the noise 
properties of catalytic-sRNA mediated regulation. 7° 

In previous studies, it was believed that the transcription 
of RNA is a standard Poisson process, which means that the 
Fano factor F, = F, = 1.8839] However, the sRNA-mediated 
gene silencing is not a pure Poisson process. Only when 
the parameter Œm is far away from the regulatory threshold, 
the Fano factor is close to 1, while the Fano factor around 
the threshold is much larger than 1 (Fig. 3(b)). The burst 
of the noise of target mRNA around the threshold is due to 
the anticorrelation between the number of mRNAs and that of 
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sRNAs."8! This conclusion can be fit to both Levine et al.’s 
model and the CSMGS model. 

Comparison of the noise properties of the two models 
suggested that the catalytic sRNA can more strictly repress 
the target mRNA at the threshold site of Q, ~ Qç. 20:46] Fig- 
ure 3(b) exhibits a right-shift threshold of the noise amplifica- 
tion effect for the CSMGS. 

Another important question that arises pertains to the 
noise transmitted in the process of sRNA-mediated gene reg- 
ulation, for example, whether the noise at the mRNA level 
would be transmitted to the protein level. Therefore, in 
the scenario of sRNA, the noise strength of protein expres- 
sion includes the bursting noise of the mRNA and protein 
production. 20! However, for protein-mediated regulation, the 


noise is calculated by F, ~ 1 +n. [20,32,44] 


4. Crosstalk between multiple target mRNAs 
mediated by a single SRNA 


Many trans-encoded sRNAs target several mRNAs, and a 
single mRNA can also be regulated by different sRNAs, which 
form a complex network.!°3449! In the case of multiple mR- 
NAs repressed by a single sRNA, mRNAs will crosstalk with 
each other through sRNA competition.>°! Considering some 
cases of double-target mRNAs, sRNA may regulate both mR- 
NAs in the stoichiometric mode, for example, RyhB under- 
goes codegradation with sodB and sdhCDAB mRNA;""! how- 
ever, in some other cases, SRNA can catalytically regulate one 
mRNA, while stoichiometrically regulating the other mRNA, 
for example ChiX represses chiP and chbB. 18-191 

Levine et al. extended the kinetic model (Eq. (3)) into a 
two-target mRNAs model and found that mRNAs exhibit hier- 
archical cross-talk that is determined by the binding rate k.!?41 
When the reporter mRNA binds with sRNA more strongly 
than another mRNA (kpr > kr), the reporter mRNA can ex- 
haust the sRNA pool quickly and release the repression of the 
later target mRNA, the function of which is similar to the 
“switch effect”, 2451 

The CSMGS model has been taken in discussing the 
crosstalk of two mRNAs by setting of one of the mRNAs reg- 
ulated in the catalytic mode and the other in the stoichiomet- 
ric mode.!! In the completely catalytic model, as ô = 1, the 
crosstalk is one-directional. The mRNA in the stoichiomet- 
ric mode exhausts sRNA to “switch on” the mRNA in the cat- 
alytic mode. This has been experimentally demonstrated in the 
ChiX—chiP-chb system.!'*!9! Before adequate induction, chb 
as a target mRNA is repressed by ChiX through base-pairing. 
On growth in chitobiose medium, the chb operon is induced. 
Along with the increase in chb, chb mRNA transforms from 
target mRNA completely repressed by ChiX sRNA to “trap- 
mRNA”, which causes the degradation of ChiX. When the 
amount of chb mRNA exceeds that of ChiX, signifying across 


the threshold depicted by Levine et al., the translation of chiP 
mRNA is released. {18-9741 

Based on the noise study of the single-target model, re- 
searchers have provided a new noise profile for two target mR- 
NAs competing to exhaust the same sRNA pool. Same as the 
single target model, of which the noise around the threshold 
(Œm ~ Os) tends to be amplified, the noise of the two-targets 
model also presents a series of peaks in the crossover region 
((Qm1 + Omi) ~ æ). PH This remarkable feature can be used 
in experiments to explore the crosstalk of double targets medi- 
ated by sRNA. 

Recently, there has been further novel progress in the 
multi-targets model where the target number is far greater than 
two.!47-50] Bosia et al. found that one target mRNA above 
its regulatory threshold is able to drive other common target 
crossover thresholds, and that hierarchical crosstalk among 


[47] Jost et al. pro- 


the multiple targets is much more notable. 
posed that sRNA weakly regulates many mRNAs as auxiliary 
targets, which have no effect on the phenotype of the organ- 
ism but provide a buffer to confer robustness to few “princi- 
pal” targets.>°! These developments have improved the un- 


derstanding of sRNA-mediated gene regulation. 


5. Sequence-function relation of sRNA in regu- 
lating target mRNA 


sRNA-mediated regulation is dependent on base-pairing 
with mRNA; therefore the functional efficiency of sRNA 
should be determined by the kinetic process of base-pairing. 
Some experiments highlighted the factors working on sRNA- 
mediated regulation, such as the length of the core interac- 
tion region, the secondary structure and the duplex formation 
energy.?°°7] The factors and the relevant energies could be 
used in the algorithms for searching the latent sRNA-mRNA 
system in vivo and were also modeled to artificially synthe- 
size the sRNA-mRNA system.3>°! Therefore, here we in- 
troduced some work involving quantification of the sequence- 
function relationship in sRNA-mediated regulation to discern 
the primary factors. 

Hao et al. constructed a series of mutants of RyhB and 
sodB, especially the complementary mutants in the core inter- 
action region, and suggested that fold-repression is exponen- 
tially correlated with the duplex energy AE in the case of no 
interruption in the core interaction region. They also found 
the energy related the secondary structure of the Hfq-binding 
region to be a good predictor for the efficacy of sRNA. [156] 
Further study suggested that the local secondary structure of 
sRNA may weaken the efficiency of sRNA owing to its un- 
availability in the kinetic process of base-pairing. According 
to the experimental results, an energetic term, AE open, was 
introduced to characterize the energy of opening of the lo- 
cal structure. We hypothesized that the competition between 
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AEopen and the base-pairing energy of the core interaction re- 
gion determined the parameters kon and koft in the first step of 
the base-pairing process. Currently, the work is under prepa- 
ration for submission. 

Many studies have emphasized the importance of the suc- 
cessive base-pairing of the core interaction region, but Peter- 
man et al. found that the loosened stem-loop structure harbor- 
ing the core interaction region tolerates mismatch in the core 
interaction region to maintain the function of sRNA.7! Peter- 
man et al. signified the binding energy of the core interaction 
and the stability of the first and third stem-loop structures of 
RyhB. A two-state binding model was built to fit the experi- 
mental data, which assumed that the relative probability of an 
unbound state to the bound state is determined by the energy 
difference AF between the bound and unbound states. 


1 


1+ e784E Me 


Panbound = 


In the heuristic model, the energy term can be written as 


Sheur = Bi (AGseea) + log(1 +exp(B2(AGsi1 — %))) 
+ log(1 + exp(f3(AGs_3 — @3))) +T. (8) 


The first term depicts the binding energy of the core inter- 
action region, and the second and third depict the stability of 
the stem-loop structures. The addictive parameters are com- 
bined into 67] 


6. Conclusion 


sRNA-mediated and protein-mediated regulations are im- 
portant for bacteria to adapt to the environment. The two regu- 
latory approaches have different features, and each of them has 
advantages and disadvantages, for example sRNA can rapidly 
respond to large changes in the input signal and shut down 
downstream expression, while protein-mediated regulation is 
better at transmitting small input signals for quantitative ad- 
justment. It is known that, sRNA functions in diverse modes 
of activating or repressing the translation of mRNA, but there 
has been less theoretical study on active regulation. Here we 
focused on the models of sRNA-mediated gene silencing. 

Two types of sRNA-mediated gene silencing have been 
demonstrated by mass equations: stoichiometric regulation, 
in which sRNA undergoes codegradation with mRNA, and 
catalytic regulation, in which partial sRNA codegradation of 
mRNA is recycled as free state sRNA. The catalytic model 
(CSMGS) inherits the linear-threshold property of the stoi- 
chiometric model, but presents stronger repression of the tar- 
get mRNA, not only in inhibiting the mRNA expression but 
also in restricting its noise strength. Although both the models 
we introduced concern the scenarios of degradation of mRNA 
caused by sRNA, Equation (3) can also be used to depict 


another kind of sRNA-mediated gene silencing where it re- 
presses the translation without degradation of mRNA, which 
also exhibits the linear-threshold property. 14 

When multiple mRNAs are regulated by a single sRNA, 
they may affect the translation of each other through sRNA 
competition. The crosstalk of mRNAs exhibits hierarchy prop- 
erty, which means the sRNA shows regulation prioritization 
among multiple mRNAs. This property has been observed in 
different experiments and models. 457! The theoretical model 
proposed that only strongly binding mRNA can efficiently re- 
lease the repression of the weak binding mRNA regulated by 
the same sRNA, but the reverse is invalid. Therefore, the 
strong binding mRNA can be considered as a switch for the 
weakly binding mRNA in both stoichiometric and catalytic 
modes. 

In the context of a single target, theoretical models 
demonstrated that the binding rate k of sRNA-mRNA inter- 
action has no influence on the threshold value, but just de- 
termines the stiffness of the transition from repression to 


[20.24] However, when multiple targets are expressed 


release. 
simultaneously, the relative binding rate contributes to the hi- 
erarchic crosstalk of mRNAs. In the last section of the pa- 
per, we have introduced some studies on the factors related 
to the efficacy of sRNA, most of which were believed to de- 
termine the binding rate of the kinetic process. Recently, Fei 
et al. applied single-molecular fluorescence in situ hybridiza- 
tion and single-molecular localization-based super-resolution 
microscopy to measure the numbers of sRNA, mRNA, and 
sRNA-mRNA complex, and then fitted them in the kinetic 
model to obtain kinetic parameters, which is the first time that 
the kinetic parameters in sSRNA-mediated regulation have been 
measured in vivo. The results suggested that the base-pairing 
of SgrS sRNA with ptsG mRNA is a reversible process with a 
large dissociation constant.*! This work inspired us to com- 
bine the measured kinetic parameters and secondary structures 
of sRNA to further study the kinetic process of base-pairing 
and the internal correlation, which will promote the develop- 
ment of the theoretical model. 


References 


[1] Masse E, Vanderpool C K and Gottesman S 2005 J. Bacteriol. 187 6962 

[2] Opdyke J A, Kang J G and Storz G 2004 J. Bacteriol. 186 6698 

[3] Bejerano-Sagie M and Xavier K B 2007 Curr. Opin. Microbiol. 10 189 

[4] Vogel J and Papenfort K 2006 Curr. Opin. Microbiol. 9 605 

[5] Brantl S 2009 Future Microbiol. 4 85 

[6] Storz G, Vogel J and Wassarman K M 2011 Mol. Cell 43 880 

[7] De Lay N, Schu D J and Gottesman S 2013 J. Biol. Chem. 288 7996 

[8] Geissmann T A and Touati D 2004 EMBO J. 23 396 

[9] Henderson C A, Vincent H A, Casamento A, Stone C M, Phillips J O, 

Cary P D, Sobott F, Gowers D M, Taylor J E and Callaghan A J 2013 

RNA 19 1089 

[10] Otaka H, Ishikawa H, Morita T and Aiba H 2011 Proc. Natl. Acad. Sci. 

USA 108 13059 

[11] Sauer E, Schmidt S and Weichenrieder O 2012 Proc. Natl. Acad. Sci. 
USA 109 9396 


128703-6 


Chin. Phys. B Vol. 24, No. 12 (2015) 128703 


[12] 
[13] 
[14] 


[15] 


[16] 
[17] 


[18] 
[19] 


[20] 
[21] 
[22] 


[23] 


[24] 
[25] 
[26] 
[27] 


[28] 


[29] 
[30] 
[31] 
[32] 
[33] 


[34] 


Aiba H 2007 Curr. Opin. Microbiol. 10 134 

Desnoyers G and Masse E 2012 Genes Dev. 26 726 

Prevost K, Salvail H, Desnoyers G, Jacques J F, Phaneuf E and Masse 
E 2007 Mol. Microbiol. 64 1260 

Soper T, Mandin P, Majdalani N, Gottesman S and Woodson S A 2010 
Proc. Natl. Acad. Sci. USA 107 9602 

Guantes R, Cayrol B, Busi F, et al. 2012 Mol. BioSyst. 8(6) 1707 
Lalaouna D, Simoneau-Roy M, Lafontaine D and Masse E 2013 BBA- 
Gene Regul. Mech. 1829 742 

Overgaard M, Johansen J, Moller-Jensen J and Valentin-Hansen P 2009 
Mol. Microbiol. 73 790 

Figueroa-Bossi N, Valentini M, Malleret L and Bossi L 2009 Genes 
Dev. 23 2004 

Hao Y, Xu L and Shi H 2011 J. Mol. Biol. 406 195 

Beisel C L and Storz G 2011 Mol. Cell 41 286 

Gogol E B, Rhodius V A, Papenfort K, Vogel J and Gross C A 2011 
Proc. Natl. Acad. Sci. USA 108 12875 

Plumbridge J, Bossi L, Oberto J, Wade J T and Figueroa-Bossi N 2014 
Mol. Microbiol. 92 648 

Levine E, Zhang Z, Kuhlman T and Hwa T 2007 PLoS. Biol. 5 e229 
Mehta P, Goyal S and Wingreen N S 2008 Mol. Syst. Biol. 4 221 
Thattai M and Van Oudenaarden A 2004 Genetics 167 523 

Elowitz M B, Levine A J, Siggia E D and Swain P S 2002 Science 297 
1183 

Swain P S, Elowitz M B and Siggia E D 2002 Proc. Natl. Acad. Sci. 
USA 99 12795 

Paulsson J 2004 Nature 427 415 

Paulsson J 2005 Phys. Life Rev. 2 157 

Golding I, Paulsson J, Zawilski S M and Cox E C 2005 Cell 123 1025 
Elf J and Ehrenberg M 2003 Genome Res. 13 2475 

Detwiler P B, Ramanathan S, Sengupta A and Shraiman B I 2000 Bio- 
phys. J. 79 2801 

Shimoni Y, Friedlander G, Hetzroni G, Niv G, Altuvia S, Biham O and 
Margalit H 2007 Mol. Syst. Biol. 3 138 


35] 


43] 
44] 


50] 
51] 


128703-7 


Ali Azam T, Iwata A, Nishimura A, Ueda S and Ishihama A 1999 J. 
Bacteriol. 181 6361 


Vogel J 2009 Mol. Microbiol. 73 737 


] Elowitz M B, Levine A J, Siggia E D, et al. 2002 Science 297 1183 


Thattai M and Van Oudenaarden A. 2001 Proc. Natl. Acad. Sci. USA 
98 8614 
Swain P S 2004 J. Mol. Biol. 344 965 


] Pedraza J M, Paulsson J. 2008 Science 319 339 


So L, Ghosh A, Zong C, et al. 2011 Nat. Genet. 43 554 


van Kampen N G 1981 Stochastic Processes in Physics and Chemistry 
(Amsterdam: North-Holland) p. 219 


Jia Y, Liu W, Li A, et al. 2009 Biophys. Chem. 143 60 
Li R, Xu L and Shi H 2015 J. Theor. Biol. 365 377 
Gillespie D T 1977 J. Phys. Chem. 81 2340 
Klironomos F D and Berg J 2013 Biophys. J. 104 951 


] Bosia C, Pagnani A and Zecchina R 2013 PLoS One 8 e66609 
] Levine E, Huang M, Huang Y, et al. 2008 URL matisse. ucsd. edu/ 


hwa/pub/sRNA-noise.pdf Preprint 

Gottesman S 2004 Ann. Rev. Microbiol. 58 303 

Jost D, Nowojewski A 2013 Biophys. J. 104 1773 

Noorbakhsh J, Lang A H and Mehta P 2013 PloS One 8 e72676 
Peterman N, Lavi-Itzkovitz A and Levine E 2014 Nucleic Acids Res. 
42 12177 


] Rodrigo G, Landrain T E and Jaramillo A 2012 Proc. Natl. Acad. Sci. 


USA 109 15271 
Busch A, Richter A S and Backofen R 2008 Bioinformatics 24 2849 
Richter A and Backofen R 2012 RNA Biol. 9 954 


Hao Y, Zhang Z J, Erickson D W, et al. 2011 Proc. Natl. Acad. Sci. 
USA 108 12473 


] Rice J B and Vanderpool C K 2011 Nucleic Acids Res. 39 3806 


Fei J, Singh D, Zhang Q, Park S, Balasubramanian D, Golding I, Van- 
derpool C K and Ha T 2015 Science 347 1371 


